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Abstract 



We study 2D-polaritons in an atomically thin dipole-active layer (2D-crystal) 
placed inside a parallel-plate wavegude, and investigate the possibility to ob- 
tain the localized wavegude modes associated with atomic defects. Consider- 
ing the wavegude width, I, as an adjustable parameter, we show that in the 
waveguide with I ~ 10 4 a, where a is the lattice parameter, the localized mode 
can be provided by a single impurity or a local structural defect. 
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I. INTRODUCTION 



The mechanism of the photon localization via the resonant coupling between photons 
and the local excitations inside the polariton gap was for the first time proposed in Refs. 
|U . Considering a dipole-active impurity atom in an isotropic frequency-dispersive medium, 
the authors discovered the photon-atom bound states. They showed that when the atom 
transition frequencies fall inside the polariton gap, the radiative relaxation of the bound 
states is suppressed, and the field is localized around the impurity. 

In the recent papers |||| , we considered another type of the impurity- induced polariton 
states associated with the phonon local states in polar crystals. It was shown that near 
the bottom of the polariton gap, the local states are predominated by the long- wavelength 
modes and have the macroscopic localization and coherence radiuses. We attributed these 
features to the singularity of the density of states at the gap bottom, which also causes the 
absence of the lower threshold for the impurity strength. The singularity is provided by the 
long-wavelength polaritons and is generic for any isotropic dipole-active phonon mode with 
the negative dispersion. 

The long-wavelength nature of the polariton states allowed us to analyze the crossover 
between the polariton and phonon local states within the continuum approximation. Our 
results show that the crossover takes place in a relativistically narrow interval near the bot- 
tom of the polariton gap 0. A small width of the crossover region and a strong suppression 
of the photon content of the polariton local states are caused by the fact that the typical mo- 
mentum of the modes dominating in these states, fc max ~ /3 1//2 a _1 , is much greater than the 
cross-resonance momentum, ko ~ where a is the lattice parameter, and (3 = v/c is the 

ratio between the phonon velocity and the speed of light. To eliminate this disproportion, 
one needs to lower the group velocity of electromagnetic waves in the active medium. It can 
be achieved if the medium is placed inside a waveguide. For instance, in the parallel-plate 
waveguide, the dispersion law of the propagating modes, Un(k) = C\J (7m/2l) 2 + k 2 , provides 
the reduction of their group velocity in the long-wavelength region. The phonon and photon 
velocities become comparable in the cross- resonance region if the spacing between plates, I, 
is of the order of 10 6 a. However, as long as the waveguide spectrum contains the activa- 
tionless mode (n = 0), there is a guarantee that fc max is far away from the cross- resonance 
point. 

In the present paper, we investigate the polarization waves in an atomically thin dielec- 
tric layer ( ionic 2D-crystal ) placed in a narrow parallel-plate waveguide. The srtuctural 
stability of 2D-crystal is provided by a strong confining potential, which also eliminates the 
activationless mode from the waveguide spectrum. A sub-micron crystal film grown on a 
substrate can be a physical realization of this model. The coherent interaction between the 
waveguide modes and the polarization waves gives rise to new exitations - 2D-polaritons. 
They form two polariton bands, and the position of the maximum in the lower polariton 
branch depends on the waveguide width. We show that A; max ~ k for / ~ /J -2 / 3 a, and fc max 
tends to (3 2 ^ a~ l when / ^> (3~ 2 ^a. In the latter case, the frequency region of the polariton 
local states is enlarged by the factor of (3~ 1 ^ ~ 10 2 , comparing to the 3D-case. The electric 
field in these states is always delocalized across the 2D-crystal within the distance of the 
order of I from the defect due to the contribution from the upper polariton band. Upon 
increase of the impurity strength, the localization length decreases and contributions to the 



2 



field from the lower and upper branches begin to compete everywhere in the waveguide. 
When the localization length becomes comparable with the waveguide width, the field is 
no longer confined near the dielectric layer and the polariton local state transforms into the 
localized waveguide mode. We estimate the corresponding value of the impurity strength, 
and show that, in the waveguide with / ~ (3~ 2 ^ 3 a, the localized mode can be provided by a 
single impurity or a structural defect. 



II. POLARITONS IN THE IONIC 2D-CRYSTAL 



Let us consider an atomically thin dielectric layer between two perfectly conducting 
sheets [Fig.l]. The layer presents a regular ionic 2D-crystal, which is stabilized in z = 
plane by a strong restraining potential. This increases the activation energy of the off- 
layer lattice vibrations and shifts the corresponding phonon modes much higher than the 
modes with the in-plane polarization. Assuming an infinitely strong confining potential we 
eliminate the off-plane phonons from our model. Among the remaining in-plane phonons, 
for the sake of simplicity, we consider only a single transverse dipole-active mode with the 
isotropic spectrum and negative dispersion. In the long-wavelength region, we can use the 
standard approximation, Q 2 (k) ~ — v 2 k 2 , where v sets the range of the typical phonon 
velocities. The dipole-active excitations in a thin layer have a high decay rate, unless induced 
electromagnetic radiation from the layer is compensated. In our model it is provided by the 
coherent coupling between the lattice excitations and the eigen modes of the waveguide. In 
the parallel plate waveguide, there are two types of propagating modes [Fig.l]. TM- modes 
(transverse-magnetic) involve electric field directed across the dielectric layer and, therefore, 
they cannot be activated in our model. Electric field in TE-modes (transverse-electric) 
is directed along the layer and, therefore, these modes can be excited along with the the 
transverse-optical phonons. Introducing the surface polarization Q(x,y) associated with 
these 2D-phonons and considering its dynamics together with the corresponding TE-modes 
of the waveguide, we obtain the following system of equations: 

c2 ^rr + - c " k2 ) E * = -4^ 2 Q k 5 (*), (i) 



ad 2 

[ L o 2 -n 2 (k)}Q k = -^-E ]l (0), (2) 

where Qk arid E^(z) are the 2D-Fourier amplitudes of the polarization and electric fields, k 
is a 2D-wave vector, d is the phonon-photon coupling parameter ( ion plasma frequency). For 
the " order of magnitude" estimates we assume in this paper that d ~ Q ~ v/a. For later 
convenience we introduce the dimensionless variables Ik — > k, Iuj / c — > uj, and IQ/c — > Q, so 
that il 2 (k) Qq — (3 2 k 2 , and f2 ~ A; ~ r]f3, where rj = l/a is a large parameter. 

The electric field in Eq.(l) is confined between the conducting plates and its normal 
derivative at the dielectric sheet has a discontinuity caused by the surface polarization. 
Solving Eqs.(l, 2) under these conditions, we obtain the following dispersion equation: 

Q 2 - uj 2 c tan Vuj 2 - k 2 

o =0 . (3) 

uj 2 V^ 2 - k 2 V ; 
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where 5 = aid 2 / 2 c 2 ~ 7](3 2 is a small parameter. 

This equation defines a series of polariton branches, oo = u n (k) , with their activation 
frequencies given by the equation [Fig. 2,3]: 

u 2 — = —5u n tanu; n . (4) 

Analysis of Eq. (3) shows that the upper (n > 1) branches have the quadratic "large mo- 
menta" asymptotes: 

a£(*)«a£ + * 2 , (5) 

where u n xs ir (n — 1 / 2) 3> Oo- Therefore, their spectral bands all overlap and form a 
common upper polariton band. 

The lower (n = 0) polariton branch is separated from the others by the polariton gap, 
the bottom of which coincides with the maximum in this branch. Evaluation of the group 
velocity at the center of the Brillouin zone, 



du 2 (k) 
dk 2 



-{? + (-^- - lV -{? + ^ , (6) 



fc=0 2u; 2 \sin2a;o 



shows that the lower branch has a negative dispersion if / < f3~ 2 ^ 3 a ~ 10 4 a. In such nar- 
row waveguides the photon-phonon coupling is negligible and the phonon branch, with its 
maximum located at k — 0, remains practically unaffected by the field. 

As / increases, the maximum moves away from k — 0, reaching the cross-resonance region, 
k ~ k , at I ~ (3~ 2 / 3 a. We restrict our further consideration to the case, / ^> /3~ 2//3 a, only. 
It guarantees that the maximum of the lower branch is located far away from k , in the 
region where k ^> Q,u>, and Eq.(3) can be approximated as follows: 

n 2 - lu 2 « -^cu 2 (cu 2 + 2A; 2 ) . (7) 

A positively defined solution of Eq.(7) gives the "large momenta" asymptote of the lower 
polariton branch, 

u^k)^^(l-y). (8) 



The corresponding dispersion curve reaches its maximum at the point 
where it sets the bottom of the polariton gap: 

-L, « fiS - 3^ (^) 2/3 ~ (il - sflv ( ^„) 2/3 £ . (io) 

In the immediate vicinity of /c max the polariton dispersion law can be presented as follows: 
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<"o(*) « o£« - 3/3 2 (k - k max ) 2 , (11) 
and the asymptote of the density of states near the gap bottom has the form: 

P(" ) °c - - 7= = = - . = (12) 



Restoring the true dimensionality of the variables and comparing our results with those 
obtained for 3D-polaritons: 



p{uj 2 ) oc 



a(a/c max ) 2 (3 

2 ~~^2 ~ Q_ . /TV Z7"^' 



one can see that the maximum of the polariton curve is now located closer to the cross- 
resonance point, ko ~ /9a" 1 , and the bottom of the gap is shifted toward the phonon acti- 
vation frequency. Also, the singularity of the density of states is strengthened by the factor 
of /3 -1 / 3 ~ 10 2 . Setting p(co 2 ) ~ 1/A, where A is the width of the polariton band, one can 
estimate the frequency range where the singularity prevails: 

A 



V<^~P 2/3J =-~P 2/3 V . (13) 

Comparing to the 3D-case, it is enlarged by the factor of /5" 1 / 3 ~ 10 2 , what broadens the 
region of the local states dominated by long-wavelength modes. 



III. DEFECT-INDUCED LOCAL STATES 

If a point-like defect is embedded in the dielectric layer, it modifies Eq.(2) as follows: 

r 2 <^2n\i^ ad 2 aa 2 Q(0) . . 

\uj 2 - n 2 (k) ] Q k = -— £ k (0)+ ^ J , (14) 

where Q(0) is the polarization of a defect and S is the total area of the layer. The strength 
of the defect, a, depends on its mobility and binding energy in a crystal. In the case of 
an isotope impurity a = —u 2 5m/m; for a non-isotope impurity or a structural defect we 
assume that a ~ a^ 2 ,, where a is a numerical parameter. 

Solving Eqs.(l,14) we obtain ( in the dimensionless variables u, A;, and Q): 

_, . . aa 2 lQ(0) uj 2 sm(Z ( 2 ^ 2 .u^tanCA" 1 . . 



where we denote Z = Z" 1 (/ ± z) and ( = \fuj 2 — k 2 . 

Using this equation one can express Qk via the polarization of the defect, Q(0), and 
obtain then the spectral equation for the local state: 



When the frequency approaches u; max , the integral diverges at the "surface" k 2 = k^ ax . Near 
w max we can approximate Eq.(16) as follows: 

oo 

1 ~ _?L (-Y x f ^ (i 7 ) 

~ 2c 2 V2tt/ 7 c 2 - cu 2 + 3/? 2 ffc - fc max ) 2 ' 







Retaining here only the singular part of the integral, we finally obtain: 

y/u? - < ax « aQ (Qfcma ^ ~ a^ 5/3 - (18) 
V max 8tic 2 PV3 

Comparing to the 3D-case (^\Joo 2 — c^ 2 iax ~ arjj3 2 ^j , the separation of the local state from 

the bottom of the gap is enlarged by the factor of (3~ 1 ^. 

Equation (18) defines the eigen frequency of the local state near the bottom of the gap. 
To evaluate the localization radius of this state one needs to consider the spatial distribution 
of the electric and the polarization fields. The inverse Fourier transformation of Eq. (15) gives 
us: 

oo 

au 2 a 2 Q(0) f dkk H (Ilk) sin (Z() ( 
[Z,T) 4tt/c 2 X y (cu 2 -ft 2 )CcosC + ^ 2 sinC ' [ ' 

— oo 

where H (Tlk) is the Hankel function of first kind, H = l~ l r . 

It follows from the spectral properties of our system, that the denominator of the inte- 
grand, considered as a function of u 2 , has a series of isolated simple zeroes: 

F (uj, k) = [(uj 2 -n 2 )( cos ( + 5uj 2 sin (] ex Y[[u 2 -ujl(k)] , (20) 

n 

where the index n enumerates different polariton branches, u 2 (k) . Since the frequency of 
the local state, u, lies in the polatiton gap, all /c-zeroes of F (u, k) are removed from the 
Re/c-axis. Closing the integration contour in Eq.(19) through the upper half-plane, we can 
calculate the integral by the method of residuals: 

E(z, r) = ^y (0) x J2 {H (TZk) sin (ZQ [F (u>, k)]~ l } kn , (21) 

where k n is a pole associated with n-th branch. 

Near the bottom of the gap (uj > cj max ) the factor presenting the lower polariton branch, 
uj 2 — uj 2 (k) , is small for k ~ k max . It suggests that the poles of [F (oj, k)]~ x corresponding 
to the lower branch are located in the vicinity of A; max . Using Eq.(ll) one can find: 



/ id 2 UJ 2 

k() — A^max i iy = ^ max (^2) 

Taking into account that /c max ~ r]f3 2 / 3 ^> co max ~ r)/3, and x ~ a /c max <C /c max for a weak 
defect, one can obtain: 
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Res [F (u,k)]~l 



dF[u (k),k] 
dk 



exp ( k maxJ (23) 



max 

+ " 3(3 2 x k 



To evaluate zeroes of F (uj, k) corresponding to the upper branches, we first impose the 
upper limitation on the width of the waveguide, such that (3~ 2 ^a <C I <C (3~ l a. In this case, 
the phonon band lies well below the waveguide cut-off frequency and all high-order polariton 
branches are practically indistinguishable from the parent TE-modes. It gives us: 

k± = ± i y/u% - uj 2 w± iuj n , (24) 

and the corresponding residuals: 

Res[F(^)]7^^. (25) 

Finally, using Eqs. (21-26), we can find the large-distance asymptote of the electric field: 



E(z,r) cx TZ- 1/2 



exp (—Tlx — zk maiX ) ^ exp (— TZu n ) sin (Zuj n ) cosc<j„ 
6 (5 2 k kUL , 5 uj 2 ool/ 2 



(26) 



This result presents the field as a sum of contributions from all branches. However, 
because u n lie close to 7r (n — |) for the upper branches, the only several terms corresponding 
to the low-located branches can be retained in the sum. 

Equation (26) explicitely demonstrates the localization of the field in radial directions, 
whereas, the ^-profile of the field depends on the defect strength, a. For a weak defect, when 
x <C ui( localization length is much greater than /), the first term of Eq. (26) dominates 
at TZ ^> x _1 and confines the field within the fc~a X -wide layer around the dielectric sheet. 
However, in the kernel of the local state, within the distance TZ < x _1 from the defect, 
next terms of Eq.(26) begin to compete with the lower band contribution. This leads to 
"derealization" of the field in z-direction within the kernel . 

Upon increase of the defect strength, the local state moves away from u max and the 
localization length, x -1 , decreases. When we are reaching a CT ~ uJi/k m3uX , the first term of 
Eq.(25) is no longer dominating in E(z, r). In this case, despite the localization of the field 
within the /-range in radial directions, its z -confinement completely disappears. Such a 
state can be qualified as a local waveguide mode. Analysis of the structure and properties of 
these states requires one to consider waveguides with / ~ (3~ 2 ^ 3 a, where k 

max ~ k{). Recalling 

the result concerning the energy distribution for 3D local polaritons 0, W /7 fieid/H / me ch ~ 
(^oAmax) 4 , we can expect the energy equipartition in the local waveguide modes. Also, 
since A; max ~ lo\ in the sufficiently narrow waveguides, there the local modes can be provided 
by defects with a < 1, such as isotop impurities or local structural defects. A more detailed 
analysis of this case will be presented elsewhere. 
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